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The antiferromagnetic Heisenberg model on the two-leg ladder with exchange interactions along 
the chains, rungs, and diagonals is studied using the Jordan- Wigner transformation and bond-mean- 
field theory. The inclusion of all three couplings introduces frustration to the system and depending 
on their relative strengths the ladder can adopt one of three possible magnetically-disordered gapped 
states. The phase diagram found in this mean-field approach is in very good agreement with the 
one calculated by Weihong and colleagues using the Lanczos exact diagonalization method. By 
analyzing the ground-state energy we study quantum criticahty when the coupling parameters are 
varied at zero temperature. We study the effect of temperature on the phase boundaries, and 
find that the system shows thermally-induced criticality for some values of the rung and diagonal 
couphng constants. All the phase transitions encountered in this system occur between disordered 
phases, and are all caused by frustration. 



I. INTRODUCTION 

Quantum criticality is a current highly debated is- 
sue in strongly correlated electron systems.— There are 
three kinds of strongly correlated systems. The first ones 
are those with only localized electrons like the quantum 
Heisenberg-type spin systems, where only the spin de- 
grees of freedom contribute to the physical properties. 
The second ones are those systems where electrons are 
mobile, and both the spin and electron degrees of free- 
dom are relevant like in the high-Tc materials away from 
half-filling. The third type of systems are those fermionic 
systems with both localized and itinerant electrons like in 
the Kondo-type (heavy- fermion) systems. In this work, 
we develop an analytical approach to study the quan- 
tum criticality phenomenon in the frustrated antiferro- 
magnetic (AF) two- leg Heisenberg ladder, and the ef- 
fect of temperature on this criticality. Because, absolute 
zero temperature cannot be reached in any experiment, 
it is important to investigate such temperature depen- 
dence. In the absence of frustration, the two-leg lad- 
der has been analyzed extensively both numerically and 
analyticallyi^ The case with second-nearest-neighbor in- 
teractions along the diagonals is of interest because this 
interaction adds frustration to the system, and there 
is a possibility that in real two-leg ladder materials it 
might be significant. Some examples of real two-leg lad- 
der materials are SrCu203f^ Cu2(C5Hi2N2)2Cl4r'i^ and 
La6Ca8Cu2404ii^ii Existing numerical data on the frus- 
trated ladder indicate that when the diagonal interaction 
is varied the system undergoes a quantum phase transi- 
tion. The main motivation for the present work is reach- 
ing better understanding of quantum criticality in the 
frustrated two-leg ladder from a microscopic point of view 
and study the effect of temperature on it. the method we 
develop is an analytical approach based on the Jordan- 
Wigner (JW) transformation. It is first tested at zero 
temperature by making sure it reproduces all the numer- 
ically derived exact existing results, then it is applied at 
finite temperature. 
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FIG. 1: The two- leg ladder showing the couplings along the 
chains, rungs, and diagonals is displayed. 



The Hamiltonian for the spin-i two-leg ladder with 
diagonal interactions is written as 
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H = J Sij • Si+i J + J_L Si,i • Si 

i j=l i 
N 

+ Jx ■ Si+i_2 + Si+i^i • Si_2), (1) 



where J is the coupling along the chains, J± the trans- 
verse coupling, and Jx the coupling along the diagonals 
as seen in Fig. [T] The index i labels the position of the 
spins along the two chains, each of which has N sites. The 
first term sums the interactions of nearest-neighboring 
spins along the chains (legs) of the ladder, the second 
term sums the interactions of the spins along the rungs, 
and the third term sums the interactions along the diag- 
onals. As usual, Sij is the spin operator. 

The frustrated two-leg ladder has been studied numer- 
ically using Ising and dimer expansions;^ Lanczos diago- 
nalization technique ^iSiiSiiS and density-matrix renormal- 
ization group (DMRG) J^iiiii^iHiii It has also been stud- 
ied analytically and we will next summarize briefly the 
different analytical theories that have been applied at 
zero temperature. The (dimerized) valence-bond spin- 
wave theory by Xiar>i^ was used to study the system. 
This theory breaks down for certain values of the cou- 
plings. For example it works for J/J± < 1/2 (so in the 
strong coupling limit with J± ^ J) when Jx/Ji_ = 0. 
The abelian bosonization technique applied by Weighong 
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et aZ4, and the non-Abelian bosonization method apphed 
by Allen and coworkersi^ work only in the the weak cou- 
pling limit Jx,J± < J- The Lieb-Mattis theorem applied 
by Hakobyar>ii puts limits on the transition line such as 
J_L < 2Jx but cannot determine its exact position. The 
reformulated weak-coupling field theory of Starykh and 
Balentsi^ works in the limit Jx,J± < J only. These 
authors found that in the quantum model the classical 
transition at J± = 2Jx splits into two, implying the 
occurrance of a new phase, which was later disproved 
numericallyFi^ The non-linear Sigma model was consid- 
ered by Nedelcu et alJ^ who focused on the case where 
the diagonal couplings are different. However, for the 
case with equal diagonal interactions there are portions 
of the phase diagram where their theory fails. The non- 
perturbative effective field theory of Cabra et al^ ex- 
amined which spin bonds are the strongest in order to 
put limits on weak and strong coupling regimes. How- 
ever, no zero-temperature phase diagram was produced 
in their work. 
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FIG. 2: In the left panel, the three possible ground states 
of the system in the Ising limit, namely the Neel state, the 
ferromagnetic chain state, and the ferromagnetic rung state 
are drawn. In the right panel, the labeling of sublattices cor- 
responding to the short-range spin orders that replace the 
long-range ones are shown for the Heisenberg limit. 



II. METHOD 



None of these analytical methods was able to create 
a complete phase diagram including the weak, inter- 
mediate, and strong coupling regimes which compares 
well with numerical data because of their various lim- 
itations. The exact numerical methods provide reli- 
able information about the states of the system, but 
in order to understand even these numerical results one 
needs to develop analytical approaches, which are read- 
ily generalizable to finite temperature, and which ap- 
ply to all coupling regimes. Using the bond- mean-field 
theory2ii^i2^i^ (BMFT) we seek better understanding of 
the phases of the frustrated two-leg ladder at both zero 
and nonzero temperatures. The BMFT is a mean-field 
theory that is based on the spin bond parameters. These 
parameters are not related in any way to any kind of 
long-range order. They are rather related to the spin- 
spin correlation function {S~ Sj'), with i and j labelling 
adjacent sites in the direction where this correlation func- 
tion is calculated. All quantities (S'f ), with a = x, y, z, 
arc zero in BMFT, implying the absence of any sort of 
long-range magnetic order. 



This paper is organized as follows. In Sec. |TT1 we 
explain how the BMFT, which is based on the two- 
dimensional (2D) JW transformation, is applied to the 
Hamiltonian for the two-leg ladder with interactions 
along the chains, rungs, and diagonals. The Hamilto- 
nian is handled and decoupled similarly to that of Refill. 
Quantum criticality, energy spectra, mean-field parame- 
ters, free energy, entropy, and specific heat are calculated 
in Sec. IIIII Also, the zero and non-zero temperature 
phase diagrams are calculated. Criticality at non-zero 
temperature is examined in Sec. IIVI Discussion of our 
results is made in Sec. |Vl In Sec. IVII conclusions are 
drawn. 



In the Ising limit, where quantum spin fluctuations 
arc absent, the magnetically-ordered ground states are 
shown in the left panel of Fig. [2l Depending on the rel- 
ative strength of the couplings, the system can be found 
in one of the three different ordered states displayed on 
the left of this figure. When Jx <^ J±, the system 
adopts the Neel state with ferromagnetic spin arrange- 
ments along the diagonals. When Jx ^ J± the system 
adopts the ferromagnetic rung state. In this case, the 
AF arrangement shifts to the diagonals, and spins on the 
rungs are forced to adopt a ferromagnetic arrangement. 
When Jx ^ J and J± » J, the system adopts a ferro- 
magnetic chain state, where the spins on the chains are 
ordered ferromagnetically, and antiferromagnetically on 
the rungs. The phase diagram in this classical limit will 
also be given later on for comparison. 

In the Heisenberg limit with full quantum fluctuations, 
we assume that in a one-to-one correspondence each of 
the magnetic orders in Fig. [2]evolves into a state which is 
characterized by only short-range order of the same kind 
as in the ordered state. So, we replace the rigid up and 
down orientations of the spins by sublattice labels A and 
B, respectively, as shown in the right panel of Fig. [2l Ac- 
cording to Ref J^, the boundaries in the phase diagram of 
the Ising limit move as a consequence of quantum fluctu- 
ations, but the number of phases stays the same, namely 
three. For this reason, we continue to label the phases 
using the same terms as in the Ising limit except that we 
now add the word type to indicate that the phases are 
not magnetically ordered but arc rather characterized by 
short-range order only. So, the phases are now called 
Neel-type (N-type) state, ferromagnetic-type rung (R- 
type) state, and ferromagnetic-type chain (F-type) state. 
For the N-type state in the absence of frustration, the 
ground-state of the two-leg ladder is a gapped disordered 
spin liquid dominated by short-range AF correlations for 
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any nonzero J_]_^ The gap increases linearly with rung 
coupling, except near zero coupling where a slight down- 
ward curvature is present. The ferromagnetic- type rung 
state is also referred to as the Haldane-phase or rung 
triplet phase because in this case the system presents 
some resemblance to the 5 = 1 chain. Precisely, in the 
case J_L = and Jx = J, the Hamiltonian becomes 
H = + S,,2) ■ (S,+i,i 4- Sj+1,2), with low-lying 

excitations identical to those of the S = 1 Haldane chain 
due to the fact that the interaction is now between the re- 
sultant of the two spins on the rungsiii The N and F-type 
states are also known as the dimerized or singlet states 
because of the formation of singlets along the rungs in 
the quantum case. The two phases belong to the same 
universality class so they are only topologically distinct 
due to the symmetry of the systemj" -! '-- - ' For the N-type 
state, in the limit = cxd and Jx = 0, the ground state 
is made of independent rung singlets. 

We use the spin arrangements on the right panel in 
Fig.[2]as a starting point. We assume that in the Heisen- 
berg limit short-range AF correlations are important 
enough to justify its use. Note that these spin arrange- 
ments are not static in the quantum limit, but fluctu- 
ate while the relative average orientations of adjacent 
spins remain the same. We will confirm that the frus- 
trated two-leg ladder is characterized by quantum criti- 
cality (criticality at zero temperature induced solely by 
quantum fluctuations), and show that it is also charac- 
terized by criticality at nonzero temperature. The finite- 
temperature phase transitions we propose for this prac- 
tically one-dimensional system are not from order to dis- 
order, but they are from a disordered phase to another 
disordered one; i.e., between states that differ by their 
short-range magnetic order only. 



The JW transformation and BMFT 



The JW transformation for the two-leg Heisenberg lad- 
der is defined as^i 



Hamiltonian ([T|) we get 



i,j ~ *J ' 
i-1 2 

4>is = ^EX!"'''^] for chain 1, 

d=0/=l 
i-l 2 

4>i,2 = T^i^^^iT-dj + riis] for chain 2. (2) 

d=0/=l 

Here i and j are the coordinates along the chain and 
rung directions, respectively. The phases (fiij are chosen 
so that all the spin commutation relations are preserved. 
The c] J operator creates a spinless fermion at site 
while annihilates one, and riij is the occupation num- 
ber operator. 

After applying the JW transformation ^ to the 
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^ 1 1 

i 

+ Jx y][{4,l(^iS - 77)(cl+l,2C«+l,2 - 7:) 



+ (cl+i_iC.»+i,i - ^)(42C»,2 - ^) 



(3) 



where 
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HxY = ^ ^[4ie^-"-=c.+i,i + 42e^-"'+^'ic.+i,2 + H.c] 



N 



4f E[4ic.2+H.c.] 



Jx 



N 



is the XY Hamiltonian of the frustrated two-leg lad- 
der. In BMFT, the interacting terms of the JW fermions 
are decoupled using the spin bond parameters. This ap- 
proximation neglects fluctuations around the mean field 
points; (O- (0))(0' - (O')) « 0, where O and O' are any 
operators which are quadratic in and This yields 



00' « (0)0' + 0(0') - {0){0'). 



(5) 



To apply BMFT we introduce three mean-field bond pa- 
rameters; Q in the longitudinal direction, P in the trans- 
verse direction, and P' along the diagonal. These can 
be interpreted as effective hopping energies for the JW 
fermions^^ in the longitudinal, transverse and diagonal 
directions, respectively: 



P = {Cijclj_^_^), 

P' = {ct+ijclj+i)- 



(6) 



Keeping in mind that there is no long-range order— so 
that (5fj ) = {cl iCi,i) — 1/2 = 0, the Ising quartic terms 
in equation ([3]) can be decoupled and simplified using the 
Hartree-Fock approximation ([5]), and the bond parame- 
ters as follows: 

+ (4iCj-(-i,i)cmc1+i,i 
-(4iCi+i,i)(c^,icJ_^i 1) 



(7) 
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FIG. 3: The lattice is subdivided into two sublattices in the 
case of the N-type state. 



FIG. 4; The phase tt ahernates along the chains and is zero 
everywhere else. The flux per plaquette is tt on average.— 



Note that when decoupled in the magnetization channel, 
the quartic terms give (^c\ -^Ci^i ~ ^ [c\^i,iCi+i.i - 5) ~ 
0, which is a consequence of the absence of magnetic long- 
range order. Then, the Hamiltonian becomes 
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i j=l 



N 



N 



i 



(8) 



Next, we write this Hamiltonian using the three different 
spin configurations in the right panel of Fig. [21 These 
configurations are instantaneous (not static) configura- 
tions in which adjacent spins in any direction keep on av- 
erage the same relative orientations with respect to each 
other, but fluctuate globally on a time scale determined 
by the strongest coupling constant so that any kind of 
long-range magnetic order is absent. These fluctuations 
are a consequence of the quantum fluctuations. The three 
competing configurations of Fig. [5] lead to three different 
quantum gapped spin liquid states, each characterized by 
its own short-range spin correlations and symmetry. 



B. Neel-type State 

In the N-type state, the spin arrangement at any time 
consists on average of antiparallel spins on adjacent sites 
in both the longitudinal and rung directions; the spins on 
the diagonals thus prefer to align parallel to each other. 
In the limit Jx <^ J, J±; spins tend to form loose spin 
singlets on adjacent sites along both the chains and rungs 
if J_L J but strong spin singlets on the rungs if J± ^ 
J ^ Jx- The word Neel is not used to indicate long- 
range order but only the fact that the spins orientations 
show short AF order. For this reason, we divide the 
lattice into two sublattices; there are cf^ fermions on 
sublattice A, and j, on sublattice B, where is 
any adjacent site to Fig. [31 

Following Refi^, we set the average phase per plaque- 
tte to be TT. We choose the configuration seen in figure [H 



which was suggested by Azzouz et al. in Refi^. In the 
latter, this configuration is used to get rid of the phase 
terms in the XY Hamiltonian Hxy\ the only effect is 
that the sign of the hopping term in the JW XY term 
becomes alternated along the chains. For the Ising term, 
we set Qij = Qe^'f''-^ where Q is site independent .-2- Here 
is the phase of the bond corresponding to Fig.[4lsuch 
that = TT or along the chains, which means that the 
Q terms alternate sign along the chains just like the XY 
terms do. This is necessary in order to recover the proper 
result in the limit Jx and J± go to 0, in which we get a 
result comparable to that of des Cloiseaux and Pearson^^ 
for the spin excitation spectrum for a single Heisenberg 
chain, namely E{k) = ^J|sinfc|. The alternating sign 
for Q can be justified in the same way as in Refi^. 

A ID Fourier transform along the chains is performed 
on the Hamiltonian keeping the chain index in the real 
space. The mean-field Hamiltonian can be expressed in 
the form 



H 



2NJQ'^ + NJ^P'^ + 2NJxP'^ (9) 



with the Nambu spinor defined by 

* — \Hk Hk ^2k ^2k J' 

and the Hamiltonian density given by 



(10) 



V. = 



with 



/ iJiSink Jxicosk 

-iJisink ^ 

Jxicosfc ^ 

\ JxiCOsfc — iJisinfc 



Ji = J(1 + 2Q), 
^±1 = J±(1 + 2P), 
Jxi = Jx(l + 2P')- 



Jx 1 COS k 
iJi sin k 

/ 



(11) 



(12) 



Diagonalizing Ti, yields the energy eigenvalues ±Eni and 
±En2 where 



ENi{k) ^ Jxicosk + \j Jl SVC? k + 



EN2{k) = Jxicos/c - \ Jls\T?k+^^. (13) 



Note that the subscript TV in E^p is used to remind our- 
selves of the N-type state. Similarly, the eigenenergies 
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of each of the remaining states will be labeled using its 
appropriate subscript. The free energy corresponding to 
our Hamiltonian is 



7 1 p2 



E E E 



j,s/3£;jvp(fe)l 



(14) 



The parameters are determined by minimizing F with 
respect to Q, P, and P', a calculation which leads to the 
following set of self-consistent equations 



Q 



EE 



8NJ 



k p=l,2 



EE 



k p=l,2 



P' = 



4iVJ_L 



dENpjk) 
dQ 

dENpjk) 
dP 

dENp{k) 



tanh 



/3ENpik) 



dP' 



tanh 



tanh 




The partial derivatives of the energies with respect to Q, 
P, and P' are given by 



dE 



Np 



dQ 

dEjsip 
dP 

QEnp 
dP' 



(-l)P+i2JJisin2fc 



J2sin2fc+ (^) 
(-l)P+ViJxi 



= 2Jx cos A:, 



with p=: 1,2. (16) 



Next, we will analyze the F-type chain state. 



C. Ferromagnetic- type chain state 

In this state the instantaneous spin arrangement is AF 
along the diagonals and rungs, but ferromagnetic along 
the chains. The key thing to realize is the fact that the 
Hamiltonian is symmetric with respect to exchanging two 
spins along the rungs at even sites^ii^ So we will get sim- 
ilar spectra if the diagonal terms alternate sign and the 
terms along the chains have all the same sign (phase tt per 
plaquette). The spectra are the same as for the N-type 
state except that Ji and Jxi are exchanged. Figure [S] 
illustrates this correspondence. 

A similar calculation was done by Xi Dai and Zhao-Bin 
Su29 for the two-leg ladder without diagonal interactions. 
They chose to alternate the index of the legs along the 
chain. They argued that this is justifiable because it 
still preserves the commutation relation when they use 
the JW transformation and because it gives them the 
expected results. If you twist their ladder to place the 
sites that are labeled 1 and 2 so that they each form a leg 
of the ladder you will realize that instead of calculating 
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FIG. 5: The coupling along the diagonal (chain) in the N-type 
state is equivalent to the coupling along the chain (diagonal) 
in the F-type state. 
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FIG. 6: The lattice is drawn in the case of the F-type state. 



the interaction along the legs of the ladder in the N- 
type state they actually consider the interaction along 
the diagonal for the F-type state. 

In their mean-field approach they end up with better 
results at J_l/J = 1 (and Jx/J = 0) than with the ap- 
proach of Azzouz et al^ which we are using here. We are 
however not using their approach because it does not de- 
scribe properly the limit J±/ J 0. On the contrary, the 
method we use here describes well this limit. The right 
description of this limit is crucial for any investigation of 
the phase diagram. 

The bipartite lattice for this state can be seen in Fig.[Sl 
We define two new chains 1' and 2' which are obtained 
by relabeling the sites along the diagonals as indicated 
in Fig. [6] In this way the Hamiltonian density we get for 
the F-type state has the same expression as the one of 
the N-type state PD) with Ji replaced by Jxi and vice 
versa. Explicitly, the energy eigenvalues are now ±Epi 
and ±Ep2, with 



J2 

Epiik) ^ Jicosk + \j J'^j^ sin^ k H — ^ , 



J2 

EF2ik) = Jicosfc- \/ Jl^sin^ k + (17) 



The equations for the free energy and the mean-field pa- 
rameters have the same forms as for the N-type state, 
(HD) and (Uni); respectively, except that now the partial 
derivatives of the energies with respect to Q, P, and P' 
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FIG. 7: The bipartite character of the lattice in the case of 
the R-type state is shown. 



are 

dEpp 
dQ 

dEpp 
dP 

dEpp 
dP' 



2 J cos k, 

(-i)f+ViJii 



Jl 1 sin^ k 



(-l)P+i2JxJxisin2fc 



(18) 



with p = 1,2. 



_L1 

4 



with the partial derivatives of the energies with respect 
to Q, P, and P' replaced by 



dE 



Rp 



dQ 

dEjip 
dP 

dEiip 
dP' 



(„l)P+i2JJisin2fc 
Jj sin^ k + J'^j^ cos^ k 



= Ji 



(-l)P+l2JxJxlCOs2fc 



Jj sin^ fc + J cos^ fc 



with p 



(22) 
1,2. 



Now that we have derived the mean-field equations for 
all three states, we solve them in order to get the 
zero-temperature and temperature-dependent phase di- 
agrams. These equations are solved numerically, except 
in the high-temperature limit where they are solved both 
analytically and numerically. Our results will be com- 
pared with existing exact numerical data. 



III. RESULTS 



D. Ferromagnetic- type Rung State 

This state is realized when the spins arrangement at 
any time is such that the spins along the diagonals are 
antiparallel and the spins along the rungs are parallel. 
The bipartite lattice for this state is displayed in Fig. [T] 
In this state both the diagonal and the chain terms link 
A and B fermions and vice versa; i.e., they are AF. Based 
on our previous arguments in Sec. IIIBl a phase of tt per 
plaquette can be used, with the mean-field phase this 
time alternating between and tt along the chains like in 
the N-type state. 

A ID Fourier transform is performed along the chains 
while keeping the chain labels in real space. The Mean- 
field Hamiltonian has the same expression as Eq. ([9]) with 
the Nambu spinor now defined by 

= (4^ 4^ 4 4) (19) 

and the Hamiltonian density given by 

/ iJisink JxiCosfcX 

Jxicoak iJisink 
— iJisinfc Jxicosfc 
\JxiCosfc — iJisinfc / 

(20) 

Diagonalizing this matrix yields the energy eigenvalues 
iEfii and ±iJij2 with 

ERi{k) = + \J Jf sin^ -|- Jxi cos^ k, 
ER2{k) = ^ - ^Jjfsm^k + Jl^cos^k. (21) 

gain, the equations for the free energy and the mean-field 
parameters have the same form as for the N-type state. 



A. Zero-temperature phase diagram 

The free (ground-state) energies of all three states are 
calculated as functions of coupling constants and com- 
pared. From thermodynamic considerations the state 
with the lowest free energy is the stable one, and when- 
ever free energies cross a phase transition takes place. 
Since only the ratios of the couplings are important we 
define ai = J±/J and a2 ~ Jx/J- The calculation was 
carried out for different sets of values of ai and a2, with 
J being the unit of energy. We found that at some val- 
ues of these couplings free energies cross, which means 
that a phase transition occurs. You can refer to Fig. [5] 
for a couple of examples. In this way we have obtained 
the phase diagram at zero temperature. The phase tran- 
sitions found here using BMFT are first-order ones for 
all values of a2- This agrees relatively well with most 
of the work done thus far by numericalSii^ and analyt- 
ical methods.— For small a2, most numerical methods 
lacked the required accuracy to determine the order of 
the transition, but DMRG calculation by Wangi^ found 
that for a2 < 0.287 the transition is of second-order char- 
acter, and for all larger values it is first order. Because 
these transitions take place at zero T as a consequence 
of varying the diagonal interaction, they can be labeled 
as quantum phase transitionsii Experimentally for a real 
material, one can vary the pressure and hope that the 
diagonal (or any other) coupling changes enough so that 
the critical region is reached. The (ai , a2)-phase diagram 
we calculated is compared to the Lanczos-technique data 
of Ref4 in Fig. (5] The agreement between the Lanc- 
zos method data and our results is very good, a fact that 
indicates that the present mean-field treatment is accept- 
able. The line at q;2 = 1 is exact and its placement is a 
consequence of the Hamiltonian symmetry with respect 
to exchanging J and Jx , which BMFT fully satisfies. 
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FIG. 10: (a) The mean-field parameters are displayed as func- 
tions of a2 for Qi = 1. There is a transition from the N-type 
state to the R-type state at 02 = 0.7. (b) The mean-field 
parameters are displayed as functions of ai for 02 = 2. There 
is a transition from the R-type state to the F-type state at 
ai = 1.5. 



FIG. 8; The zero-T free (ground- state) energies as calculated 
by BMFT for the N-type state (N), R-type state (R), and the 
F-type state (F) are plotted versus 02 = Jx/J- (a) Qi = 
J^/J = 0.5; there is a transition from the N-type state to 
the R-type state, (b) ai = 2; there is a transition from the 
N-type state to the F-type state. There is a discontinuity in 
the free energy of the R-type state due to a sudden change 
in the bond parameters. This is of no interest because the 
only stable state is the one with the lowest free energy, and 
the transition is determined by the crossing of the lower free 
energies. 
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FIG. 9: The (ai, a2)-phase diagram we calculate is shown. It 
is compared to the Lanczos-method,— and Ising Expansioi*^ 
results. The boundaries are between the N-type state (N), 
R-type state (R), and the F-type state (F). For comparison, 
the phase boundaries (dotted lines) in the Ising limit are also 
shown. 



B. Mean- field parameters 

The mean-field bond parameters Q, P, and P' are not 
order parameters in the conventional way because they 
are not related in any way to any sort of long-range order. 
In BMFT, they are a measure of the AF fluctuations in 
the system. Each of these parameters can be interpreted 
as indicating the presence of a strong spin bond in the 
spatial direction in which this parameter is found to be 



significant. The spin bond consists of a rcnormalized 
spin singlet formed on adjacent lattice sites. It is there- 
fore important to know the coupling dependence of these 
parameters. We found that it is the combination of how 
these parameters and free energies depend on coupling 
constants that determines the phase boundaries between 
the three possible states; N-type, R-type, and F-type. 
The zero-T mean-field parameters are plotted in Fig. [10] 
as functions of 012 and ai. From the analysis of free en- 
ergy we determine the initial and final states as well as 
the transition points. In each state the parameter that 
is zero corresponds to the direction with ferromagnetic 
arrangement. For example, in Fig. Iior al for small a2^ 
we are in the N-type state with P' = and other pa- 
rameters [Q and P) hardly changing as a2 increases. At 
Q!2 = 0.7, there is a transition to the R-type state, which 
is accompanied by a sharp change where P vanishes and 
P' increases sharply. 



C. Energy gap and spectra 

Another quantity of significant importance is the en- 
ergy gap, which characterizes all three states. As a con- 
sequence all spin correlations span regions of size of the 
order of the reciprocal of the gap. We calculated the 
function of a2 for ai = \ and reported it in 
Fig. [TTJ Our result is compared to the exact DMRG 
onesi^ii^ in the same figure. One can note that the gap 
we calculate is not in good quantitative agreement, but 
as far as trends are concerned good qualitative agreement 
is found. Our gap behaves linearly in all the mean-field 
calculations whereas in the DMRG it is nearly horizontal 
for small and large values of a2 ■ The most important fea- 
ture that both BMFT and DMRG results share is that 
the gap shows a minimum at the critical value of a2 ■ In 
the DMRG resulta^^, the transition was interpreted to 
be second order because the gap vanishes (within uncer- 
tainty) at the transition. In BMFT the transition is not 
gapless but the gap becomes very small. For a2 = 0.2 the 
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FIG. 11; (a) The energy gap calculated for qi = 1 as a func- 
tion of 02 is displayed. The solid line is the data calculated 
using BMFT and the dotted one is the DMRG data from 
RefM. (b) The energy gap plotted for ai = 0.5 as a function 
of a2. The solid line is the data calculated using BMFT and 
the dotted one is the DMRG data from Ref.i^. (c) The energy 
gap is plotted as a function of ai for 02 — 0.2. 



FIG. 12: (a) The energy spectra calculated within BMFT 
(dotted lines) are compared to those from the dimer expansion 
method of Ref.— . (a) ai = 1 and 02 = 0. (b) ai — 1 and 
a2 = 0.6. 




minimum gap value from DMRG^^ is 0.004J ± 0.004J, 
and in BMFT it is 0.009J. One should stress however 
that it is the crossing of the free (or ground-state) en- 
ergies that determines the transitions in BMFT not the 
vanishing of the energy gap. 

The energy spectra for ai — I and a2 ~ 0, and for 
a2 = 0.6 both in the N-type state are shown in Fig. [T^ 
For ai = 1 and zero diagonal coupling, the exact spec- 
trum calculated numerically starts off asymmetric with 
a minimum at fc = n^n^ As frustration increases the 
local minimum at decreases rapidly but the absolute 
minimum at tt decreases only slowly. So the energy gap 
decreases slowly while the minimum is at tt, but starts 
to decrease rapidly once the absolute minimum shifts to 
k = 0. 

Within the BMFT, the shape of the spectrum for 
ai = 1 and a2 = is different from the exact one. 
But as discussed in Ref<^, the important feature shared 
by both results is the presence of an energy gap. Also, 
BMFT rightfully describes the physics of the two-leg lad- 
der in this limit, namely that the ground state consists of 
the formation of renormalized spin singlets on the rungs. 



FIG. 13: Comparison of spectra for Q2=0. 6. The dashed curve 
is the result of the dimer expansion in Ref.—. 

The difference in the curves for the energy gap (Fig. [TT|) 
is due to the complex behavior of the spectra shape as 
the couplings are varied. Near the boundary line between 
R and N-type states, for ai = 1 and a2 = 0.6 the BMFT 
spectra yield a low-lying excitation spectrum that over- 
all behaves like the dimcr-expansion data of Ref4; see 
Fig.m 



D. Nonzero-temperature phase diagram 

Unlike exact diagonalization methods, the present ana- 
lytical approach can be readily used to analyze the effects 
of temperature on the system. We repeated the same ap- 
proach as in Sec. IIII Al bv comparing the free energies of 
the three phases, this time, at different temperatures for 
various sets of coupling values. We deduced the temper- 
ature dependence of the phase boundaries in the phase 
diagram. The result of such a calculation is reported in 
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Fig. [Ml We found that as temperature increases the R- 
type state decreases in size. The sizes of the N-type and 
F-type phases increase with temperature. 

In the R-type state at zero temperature, the spins on 
the rungs arrange themselves ferromagnetically; notice 
that on any rung the pair of spins fluctuate together be- 
tween the up and down spin orientations while the two 
pairs of spins on the adjacent rungs fluctuate in the op- 
posite direction. So, the system is neither ordered fer- 
romagnetically nor antiferromagnetically. As mentioned 
earlier, the parameter P is a measure of the AF cor- 
relations along the rung direction. Because of the fer- 
romagnetic orientation, their AF correlations are zero, 
leading to P = 0. This parameter becomes non-zero as 
temperature increases as seen in Fig. [TSl because thermal 
fluctuations allow the rungs to adopt sometimes the AF 
arrangement. Note that eventually, the AF correlations 
diminish in the very high temperature limit, a fact that is 
indicated by P decreasing as 1/T after reaching a maxi- 
mum. In the Ising limit, the phase boundary between the 
R-type and N-type phases is ai = 2a2, and the bound- 
ary between the R-type and and F-type phases is ai = 2. 
Quantum fluctuations cause these phase boundaries to 
move toward the R-type state. Including thermal fluctu- 
ations seems to have the same effect. It is found that as 
temperature increases the phase boundaries move toward 
ai ~ Oil and ai = 1, respectively. As temperature rises 
P approaches but never becomes larger than P' and Q 
for ai ^ a2 = 1, Fig. [151 
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FIG. 15: The mean-field parameters are plotted versus tem- 
perature (solid lines). The dotted lines are from the high- 
temperature limit equations (|23|) . (a) 02 = Q;i = 1 in the 
R-type state, (b) Q2 = 0.5 and ai = 2 in the N-type state. 



equations are obtained 



J 



E. High-temperature regime 



P 
P' 



, keT > J, 



8kBT{l 



Jx 



■, ksT > 



ikBT 



(23) 



In the limit ksT 3> J, J± and Jx the mean-field equa- 
tions can be solved analytically. The Fermi-Dirac factors 
(tanh functions) can be expanded to first order in PEp in 
the mean-field cquationsi^ The approximation is subbed 
into Eqs. ([15]) with J2k ~* / H' ^'^'^ following set of 
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FIG. 14: Surface plot of the phase diagram showing the tem- 
perature dependence. 



We find that these equations arc independent of the state 
in which they arc calculated; i.e., whether we use the set 
of equations for the N-type, R-type, or F-type state we 
always get the same result (|23l) in the high-T regime. In 
this regime, the parameters decrease following a Curie- 
Weiss T-^ law but never vanish, excluding in this way 
the occurrence of any finite-temperature phase transi- 
tion from a state with finite spin bond order to a high-T 
state with zero spin bond order. Note that Eqs. fit 
very well the numerically calculated parameters as seen 
in Fig. [151 It is interesting to note that all the parame- 
ters have the same form and show the same dependence 
on the ratio of the coupling constant, in the direction 
in which the parameter is calculated, and temperature. 
Note that the smallest parameter corresponds to the di- 
rection in which the spins are ferromagnetically arranged, 
e.g. P is the smallest parameter in the R-type state. At 
high temperature, it is easy to see from ([^S]) that the 
parameter corresponding to the smallest coupling value 
will be smallest. So, in the high-tcmpcrature limit the 
largest coupling value determines the state of the sys- 
tem. This is why the boundaries in the phase diagram 
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FIG. 16: The parameters are plotted versus temperature for 
two coupling sets, (a) 02 = 0.5 and ai = 0.6. The phase 
transition happens at ksT / J = 0.39 from the R-type to the 
N-type state. (b) a2 = 2 and ai — 1.25. The phase transition 
happens at ksT/J — 0.64 from the R-type to the F-type 
state. 
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IV. CRITICALITY AT NONZERO 
TEMPERATURE 

For sets of coupling values (ai, 012) within the shaded 
region of Fig. [TH the thermal fluctuations can cause a 
first-order phase transition from the R-typc state to the 
F-type state or N-type state. There are no thermally- 
induced transitions between the F-type and N-type states 
because the boundary between the N and F-type states is 
not temperature dependent due to the symmetry of the 
Hamiltonian with respect to exchanging J and Jx terms. 
Note that thermally-induced transitions are not from a 
disordered phase to an ordered one, or vice versa. The 
mean-field bond parameters arc displayed in Fig. [TB] as 
functions of temperature for two sets of couplings ai and 
Q!2 such that Fig. 116b shows a transition from the R-type 
to the N-type state and Fig. [TBb shows a transition from 
the R-type to the F-type state. The corresponding free 
energies are shown in Fig. \17\ entropies in Fig. [181 and 
specific heats in Fig. [12] for the same sets of couplings. 

The entropy is calculated using 

k p=l 

+{1 - nF[E.pik)]}\n{l - nplEpik)]}} (24) 

which is derived from 5 = In Eq. Ep{k) 

refers to the energy spectra of the state where S is cal- 
culated. The specific heat is calculated using C = T^. 
The entropy shows a discontinuity at Tp, implying that 
the transition is first-order in character. At very high 
temperature entropy saturates as expected to a value of 
fc_Bln2. In Fig. [TSfa) for a2 = 0.5 and ai = 0.6, the 
phase transition happens at ksT/J = 0.39 from the R- 
type to the N-type state when T increases. In Fig. flM b) 
with a2 = 2 and ai — 1.25, the phase transition happens 
at ksT/ J = 0.64 from the R-typc to the F-type state. In 



FIG. 17: The free energies are plotted as functions of temper- 
ature. The phase transitions happen where the free energies 
cross, (a) a2 = 0.5 and ai = 0.6. The phase transition 
happens at ksT/J = 0.39 from the R-type to the N-type 
state. (b) a2 = 2 and ai = 1.25. The phase transition hap- 
pens at ksT / J = 0.64 from the R-type to the F-type state. 



these figures, the dashed lines simply indicate the tran- 
sitions between the different phases corresponding to the 
set of couplings used. We found that in the limit ai and 
Q!2 — > the jump in entropy goes to zero. For small 
coupling values this jump could be smaller than exper- 
imental precision (if a real material existed) so that it 
would become difficult to assert that the transition is a 
first-order one. For example, for a2 = 0.2 and ai = 0.24 
the jump in entropy is about 0.002/cb only. Note that 
because all three states are gapped, both entropy and 
specific heat show an activated behaviour in the viscin- 
ity of zero temperature. 

All the phases of the system arc disordered; i.e., none 
of them is characterized by long-range magnetic order of 
any kind. The proposed thermally-induced criticality can 
be seen as a remnance of the zero-T (quantum) critical- 
ity because of the temperature dependence of the phase 
boundaries. As mentionned earlier in the introduction, 
the mean-field parameters Q, P, and P' do not repre- 
sent any kind of long-range order. The zero-temperature 
phase transitions we analyzed in Sec. IIII AI all occur be- 
tween disordered phases that differ only by the way the 
spins arrange themselves on very short distances (refer to 
Fig. [2]). Therefore, we believe that the finitc-T transitions 
we find here for the present practically one-dimensional 
system are not an artifact of the mean-field character of 
BMFT. These phase transitions are due to frustration; 
i.e., they disappear once frustration is brought to zero. 
They are also a consequence of the fact that the zero- 
T quantum phase transition boundary depends on tem- 
perature. Because zero temperature cannot be reached 
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FIG. 18: The entropy is plotted versus temperature. The 
entropy is discontinuous at the transition points. The dashed 
lines indicate the transition between the different phases, (a) 
02 = 0.5 and ai — 0.6. The phase transition happens at 
ksT/J — 0.39 from the R-type to the N-type state, (b) a2 = 
2 and ai = 1.25. The phase transition happens at ksT/J = 
0.64 fr'om the R-type to the F-type state. 



0.4 




FIG. 19: The specific heat is calculated by C = Tff . The 
infinite peaks at the transitions, which are a consequence of 
the discontinuity in S{T), are not displayed for clarity, (a) 
a2 = 0.5 and ai — 0.6. The phase transition happens at 
ksT/J — 0.39 from the R-type to the N-type state, (b) Q2 ~ 
2 and ai = 1.25. The phase transition happens at ksT / J = 
0.64 from the R-type to the F-type state. 



in practice, a quantum phase transition cannot in fact 
be observed directly. For the present system, the signa- 
ture for such a transition would be the observation of the 
finite-T transitions (if an experimental system existed). 



V. DISCUSSION 

While the quantum criticality we find in this system 
has been already found elsewhere using exact numeri- 
cal methods, the criticality at finite temperature we re- 
port on here remains to be confirmed by other theoreti- 
cal methods. Experimentally, if ever a frustrated two-leg 
ladder where the strength of frustration is as important 
as the coupling along the rungs and chains existed, then 
thermodynamical measurements would either confirm or 
refute our claims. Perhaps the application of a high pres- 
sure on a two-leg ladder material could increase the diag- 
onal interaction and allow the search for this thermally- 
induced criticality. In the absence of diagonal interaction 
(frustration), the AF Heiscnberg two-leg ladder shows 
neither quantum nor classical criticality. Interestingly, it 
is the frustration that is responsible for both criticalities. 
A strong argument in favor of the existence of the finite-T 
criticality is the existence of the zero-T one itself because 
in both cases each phase boundary separates two of the 
same three disordered phases. Denying the finite-T criti- 
cality would amount to denying the zcro-T one according 
to the present theory. 



VI. CONCLUSION 

In this work, we studied quantum as well as classical 
criticality in the two leg-ladder with exchange interac- 
tions along the chains, rungs, and diagonals using the 
Jordan- Wigner transformation and the bond-mean-field 
theory. The zero-temperature phase diagram of this sys- 
tem is calculated. We found it to exhibit three quan- 
tum phases, characterized all by an energy gap and ab- 
sence of magnetic order. These states are labeled Neel- 
type (N-type), Rung- type (R-type or Haldane-type) and 
Ferromagnetic-type chain (F-type) states. This result 
agrees well with existing numerical data. The transitions 
from one phase to any of the two others are all first-order. 
Because they occur at zero temperature, they enter un- 
der the category of quantum phase transitions. When 
temperature increases for some sets of coupling values, 
the system undergoes a phase transition from the R-type 
state to the N or F-type state at a finite temperature. 
The finite temperature phase diagram is calculated as 
well. In it, the size of the R-type state becomes smaller 
while the F-type state and the N-type state increase in 
size with increasing temperature. The good agreement 
between our results and existing exact results for the 
zero-T phase diagram suggests that the present mean- 
field treatment is acceptable. The various phase transi- 
tions found in this work occur between magnetically dis- 
ordered states. It is solely the frustration, the spin bond 
parameters and the nature of short-range magnetic corre- 
lations that determine the nature of all three phases that 
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characterize the frustrated antiferromagnetic Heiscnberg 
two-leg ladder. In the present mean-field-type approach, 
all phase transitions are first-order in character. While 
at zero temperature existing exact numerical data seem 
to indicate that this is the case, at finite-temperature the 
degree of the transitions remains to be confirmed. 
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